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SlM»LlFieO SULUnON FOR STRESSES ANU OEFORNAflONS 
Bern«rd J. H«mrock and David E. Brewe* 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 
SUNHARY 

A shortcut to the classical Hertzian solution for iQcal stress and de> 
formation of two elastic bodies In contact Is presented. The shortcut Is 
accomplished by using slinpllfled forms for the elllpticlty and for the com- 
plete elliptic Integrals of the first and second kinds as a function of the 
geometry. Thus the Interdependence of these variables can be uncoupled, and 
the resulting transcendental equation, which must be solved through use of 
the computer or design charts, avoided. 

Simplified formulas that make the elastic deformation at the center of 
contact easy to calculate have been previously reported by the authors. 
However, the range of applicability was limited to el Uptlcltles greater 
than or equal to 1. This paper extends the range of validity to Include 
el Uptlcltles less than 1, that Is, where the semimajor axis In the ellipti- 
cal contact lies In a direction parallel to the rolling direction rather 
than being perpenolcular as In previous studies. Furthermore an auxiliary 
shear stress parameter Is expressed In simplified form as a function of the 
geometry. This enables a shortcut calculation to be made for the location 
and magnitude of the maximum subsurface shear stress. 


•Propulsion Laboratory, AVRADCOM Research and Technology Laboratories. 



INTRODUCTION 

The classical Hertzian solution for Oefonnation requires the calcula- 
tion of the ellipticity parameter k and the complete elliptic integrals of 
the first and second kinds, ^ and /, respectively. Simplifying expres- 
sions for k, and / as a function of the radius ratio a were pre- 
sented by Brewe and Hamrock (1977) using a curve-fit analysis. With these 
expres<5ions the deformation at the center of the contact « could be deter- 
mined, wrth a slight sacrifice in accuracy, without involved mathematical 
methods or the use of design charts. The sin^lifying expressions were use- 
ful for radius ratios ranging from circular point contact to a near line 
contact normal to the rolling direction ( i .e. , 1 .0 < a < 35) . However, 
the»*e are a number of applications for which the semimajor axis in the 
elliptical contact lies in a direction parallel to the rolling direction, 
resulting in a < 1. In local deformation due to asperity-asperity interac- 
tion, the radius ratio can range from much less than 1 (Patir and Cheng, 
19/8) to infinity. For many gear and rolling-element bearing applications 
the "run-in" surface becomes anisotropic and the radius ratio is generally 
less than 1 for local asperity contact. Other exan^les where a may be 
less than 1 are (1) Navikov gear contacts, (i.) locomotive wheel-rail con- 
tact, and (3) roller-flange contact in an axially loaded roller bearing. 
Inerefore the elliptical-contact deformation and stresses to be presented in 
tnis paper are applicable for any contact ranging from something similar to 
a disk rolling on a plate (radius ratio o • 0.03) to a ball-on-plate con- 
tact (a - 1) to a contact approaching a nominal line contact (a ♦ 100) such 
as a barrel-shaped roller against a plate. 

Thus far, we have limited our consideration to the use of the simpli- 
fied formulas in determining the elliptical-contact deformation. In this 
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paper we further illustrate their applicability in calculating surface 
stress as well as subsurface stress, which is important to the determination 
of fatigue life in rolling-element bearings. 
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SYMBOLS 

solid A 
solid 6 

diameters of contact ellipse along x and y directions. 


respectively (cm) 

2 

modulus of elasticity (N/cm ) 



complete elliptic integral of second kind 
approximate val’je of / using curve-fit equation 
applied load (N) 

complete elliptic integral of first kind 
approximate value of ^ usir.g curve-fit equation 
(Dy/0^), ellipticity 

approximate value of k using curve-; it equation 
effective radius of curvature in x-plane (cm) 
effective radius of curvature in y-plane (cm) 

(1/R + 1/R^.). curvature sum (cm~^) 

X y 

principal radii of solid a (cm) 
principal radii of solid b (cm) 
auxiliary parameter 

auxiliary parameter using curve-fit equation 
location along x-axis (rolling direc^.ion) of maximum subsurface 
shear stress (cm) 

calculated value of Xq using approximate formula (cm) 
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depth of maximum subsurface shear stress In x -2 plane (cm) 

7q calculated value of using approximate formula (cm) 

r curvature difference 

6 elliptical-contact deformation at center of contact (cm) 

T elliptical-contact deformation calculated by using approximate 

formulas (cm) 

a maxir.ium Hertzian stress (N/cm^) 

maximum Hertzian stress calculated by using approximate formulas 
(N/cm^) 

Tjj maximum subsurface orthogonal shear stress (N/cm*^) 

Tq maximum subsurface orthogonal shear stress using approximate 

2 

formulas (N/cm ) 

CONFORMAL AND NONCONFORMAL SURFACES 
Hydrodynamic lubrication is generally characterized by surfaces that 
are conformal. That is, the surfaces fit snugly into each other with a high 
degree of geometrical conformity, so that the load is carried over a rela- 
tively large area. Furthermore the load-carrying surface area remains 
essentially constant while the load is incr^-ased. Fluid-film journal and 
slider bearings are conformal surfaces. i.< journal bearings the radial 
clearance between the shaft and the bearing is typically one-thousandth of 
the shaft diameter; in slider bearings the inclination of the bearing sur- 
face to the runner is typically one part in a thousand. 

Many machine elements have contacting surfaces that ^ not confonn to 
each other very well. The full burden of the load must then be carried by a 
very small contact area. In general, the contact areas between nonconformal 
surfaces enlarge considerably with increasing load but are still small com- 
pared with the contact areas between conformal surfaces. Some examples of 
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these nonconformal surfaces are mating gear teeth, cams and followers, and 
rolling-element bearings. 

The load per unit area in conformal bearings is relatively low, typi- 
cally only 1 MN/m and seldom over 7 MN/m . By contrast, the load per 
unit area in nonconformal contacts, such as those that exist in ball bear- 
ings, will generally exceed 700 MN/m even at modest applied loads. These 
high pressures result in elastic deformation of materials such that the 
elliptical contact areas are formed for load support. The present paper 
develops simple solutions for the stresses and deformations in nonconformal 
contacts. 


CURVATURE SUM AND DIFFERENCE 


The undeformed geometry of contacting solids can be represented in gen- 
eral terms by two ellipsoids. The two solids with different radii of curva- 
ture in a pair of principal planes (x and y) passing through the contact 
between the solids make contact at a single point under the condition of 
zero applied load. Such a condition is called point contact and is shown in 
figure 1, where the radii of curvature are denoted by r's. It is assumed 
throughout the paper that convex surfaces, as shown in figure 1, exhibit 
positive curvature and concave surfaces, negative curvature. Therefore, if 
the center of curvature lies within the solid, the radius of curvature is 
positive; if the center of curvature lies outside the solid, the radiu: of 
curvature is negative. It is important to note that, if coordinates x 
and y are chosen such that 


J-> JL^_L 

"ax "bx - "ay "by 


( 1 ) 


coordinate x then determines the direction of the semiminor axis of the 
contact area when a load is applied and y, the direction of the semimajor 
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axis. The direction of motion is always considered to be along the x-axis. 
For those situations in which the principal curvature planes of the two con- 
tacting bodies are not coincident, refer to Timoshenko and Goodier (1970). 

The curvature sum and difference, which are Quantities of some impor- 
tance in the analysis of contact stresses and deformation, are 


1 

IT 



1 


( 2 ) 



(3) 


where 


1 


r 


X 


- 1 . 
*^ax *"bx 


(4) 


1 „ JL ♦ -L 

^ ""ay '’by 


(5) 


Equations (4) and (5) effectively redefine the problem of two ellipsoidal 

solids approaching one another in terms of an equivalent ellipsoidal solid 

of radii R and R approaching a plane. 

X y 

The radius ratio a is defined as 

( 6 ) 



Thus, if equation (1) is satisfied, then a is greater than or equal to 1; 
and if it is not satisfied, a is less than 1. 

SURFACE STRESSES AND DEFORMATION 

When two elastic solids are brought together under a load, a contact 
area develops, the shape and size of which depend on the applied load, the 
elastic properties of the materials, and the curvatures of the surfaces. 
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Mhen the two solids shown in figure 1 have a normal load applied to them, 
the shape of tne contact area is elliptical. It has been common to refer to 
elliptical contacts as point contacts, but since this paper deals mainly 
with loaded contacts, the term elliptical contact is adopted. For the 
special case where r^^ . r^^ and r^^^ - r^^^, the resulting contact 
is a circle rather than an ellipse. Where r^^ and r^^^ are both in- 
finite, the initial line contact develops into a rectangle when load is 
applied. 

The ellipticity parameter k is defined as the elliptical-contact 
diameter in the y-direction (transverse direction) divided by the 
elliptical-contact diameter in the x-direction (direction of motion) or 

Hh/ (7) 

X 


If equation (1) is satisfied and a > 1, the orientation of the contact 
ellipse will have the major diameter transverse to the direction of motion, 
and consequently k > 1. Otherwise, the major diameter would lie along the 
direction of motion with both « < 1 and k <1. To avoid confusion, the 
commonly used solutions to the surface deformation and stresses are pre- 
sented only for the case in which a > 1. The simplified solutions are pre- 
sented and then their application for o < 1 is discussed. 

Harris (1%6) has shown that the ellipticity parameter can be written 
as a transcendental equation relating the curvature difference (eq. (3)) and 
the elliptic integrals of the first ^ and second i kinds as 


k 


2.T- /U ♦ r)^'^ 


(B) 


where 
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A one-point Iteration method that was adopted by Hamrock and Anderson (1973) 
can be used to obtain the ellipticity parameter, where 


^n+1 


a k.. 


( 11 ) 


The iteration process is normally continued until differs from 

k^ by less than 1x10'^. Note that the elliptidty parameter is a func- 
tion of the radii of curvature of the solids only: 

• fC’a«-''bx’''ay''by) 

That is, as the load increases, the semi-axes in the x and y directions 
of the contact ellipse increase proportionately to each other, so the ellip- 
ticity parameter remains constant. 

Figure 2 shows the ellipicity parameter and the elliptic integrals of 
the first and second kinds for a range of the curvature ratio 
usually encountered in concentrated contacts. 

wnen the ellipcity parameter k, the normal applied load F, Poisson's 
ratio V, and the modulus of elasticity E of the contacting solids are 
known, the major and minor axes of the contact ellipse and the maximum de- 
formation at the center of the contact can be written from the analysis of 
Hertz (1881) as 



( 14 ) 


where 



(15) 


(16) 


In these equations, 0 and D are proportional to and a is 

y ^ 

2/3 

proportional to F . 

The maximum Hertzian stress at the center of contact car also be deter- 
mined by using equations (13) and (14) and 

y ^ 

SUBSURFACE STRESSES 

Fatigue cracks usually start at a certain depth below the surface in 
planes parallel to the direction of rolling. Because of this, special 
attention must be given to the shear stress amplitude in this plane 
Furthermore a maximum shear stress is reached at a certain depth below the 
surface. The analysis used by Lundberg and Palmgrr?n (1947) is used here to 
define this stress. 

The stresses are referred to a rectangular coordinate system with its 
origin at the center of the contact, its z-axis coinciding with the interior 
normal of the body considered, its x-axis in the direction of rolling, and 
its y-axis perpendicular to the rolling direction. In the analysis that 
follows, it is assumed that y » 0. 


9 



From Lundberg and Palmgren (ly4/) the following equations can be 


written: 


zx 


» 


tan‘"v 


cos^^ sin sin t 
^ - 0^ cos^ 


0 


( 18 ) 


X 


1 

7 




tany 


sin 0 


(19) 


z 



tan Y cos t> 


( 20 ) 


where ^ and y auxiliary parameters used in place 'if the coordinate 
set (x,z). They are defined so as to satisfy the relationship for a con- 
focal ellipsoid to the pressure ellipse (for further details see Hertz, 
i88i, and Lundberg and Palmgren, 1447). The maximum shear stress amplitude 
is defined as 


T 



max 


The amplitude of the shear stress 



~W 



is obtained from 

• 0 

- 0 


For the point of maximum shear stress 

tan^P > t 
tan'^Y - t - 1 

0 2 

IT- - l)(.t - 1) 

y 


( 21 ) 

( 22 ) 

(23) 
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The position of the maximum point Is determined by 


2q - C* 2- 
Xq . * n* 2 " 


where 


(t ♦ 1) v?:i - r 


t ./TnT 

^ “ m: r ?rnr 


Furthermore the magnitude of the maximum shear stress is given by 


(24) 

(25) 


(26) 


(27) 


V2i -T 

^0 " ®max 2 1 ( t + I ) 


(28) 


It should be emphasized that Tq represents the maximum hal^-ampl Itude of 

the subsurface orthogonal shear stress and is not to be confused with the 
# 

maximum subsurface shear stress that occurs below the center of contact on 
the plane oriented 45* to the surface. The Lundberg-Palmgren prediction of 
fatigue life Is based on the calculation and was limited to cross 

sections lying in the plane of symmetry o tne roller path (y « 0). 

SIMPLIFIED SOLUTIONS FOR a > 1 

The classical Hertzian solution presented in the previous section re- 
quires the calculation of the eiHpIcity parameter k and the complete 
elliptic integrals of the first and second kinds ^ and /. This entails 
finding a solution to a transcendental equation relating k, and / to 
the geometry of the contacting solids, as expressed In equation (8). This 
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is usually accomplished by sume iterative numerical procedure, as described 
by Hamrock and Anderson or with the aid of charts, as shown by Jones 

(1946). 

Table 1 shows various values of radius-of-curvature ratios and corre- 
sponding values of k, and / obtained from ♦he numerical procedure 
given in Hamrock and Anderson (1973). For the set of pairs of data [(k^, 
a^), i - 1, 2, ...» 26j, a power fit using a linear regression by the 
method of least squares resulted in the following equation: 

IT . for a > 1 (29) 

The asymptotic behavio? of / and ^ (a ♦ 1 implies i * ^ * «/2, and 
a ^ m implies • and / ■» 1) was suggestive of the type of functional 
dependence that / and ^ might follow. As a result, an inverse and 
logarithmic curve fit was tried for J and respectively. The following 
expressions providea excellent curve fits: 

/-1 + S for a>l (30) 

where 

q - J - 1 (31) 

and 


^ ♦ q In a for a > 1 (32) 

Values of k, and JF are presented in table 1 and compared with the 
numerically determined values of k, /, and F, Table 1 also gives the 
percentage of error determined as 
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e s 


(33) 


(I - 2)100 
z 


where 


2 » {k, /. 


(34) 


z - y} 


(35) 


Table 2 shows various values of radius-of-curvature ratios and corre- 
sponding values of 0 , D , , and o obtained exactly. Also 

y X inax 

shown in table 2 are the appropriate values 0^, 0^^, ®mdx* ^ 
obtained from using equations (30) to (32) in conjunction with equations 
(13) to (15) and equation (28). The percentage of error as determined by 
equation (33) is also given in table 2. The agreement between the exact and 
approximate solutions is indeed quite good. 

Table 3 shows various values of radius-of-curvature ratios and corre- 
sponding values of the auxiliary parameter t used in calculating the posi- 
tion and value of maximum subsurface orthogonal shear stress. For the set 

of pairs of data [(t-, a-), i * 1, 2 44] the following simplified 

formula was obtained: 


t 


1 + 0.16 csch ^ 


(36) 


The position and value of maximum subsurface orthogonal shear stress corre- 
sponding to the auxiliary parameter are giver, in table 3. The percentage of 
error e is given for each of these values. The agreement between the 
exact and the approximate values of t is quite good. Once the value of 
the auxiliary parameter t is determined, the position and value of maximum 
suDsurface orthogona’ shear stress can readily be calculated. 
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SWLU-'ltl) SOLUriUNS J-'UR U.Oi <« < i 
Table 4 gives the simplified equations for conditions where 
O.Ol < a 1 100. Recall that a ^ 1 implies k ^ 1 and equation (1) is satis- 
fied and that 0 < « < 1 implies 0 < k < 1 and equation (1) is not satis- 
fied. It IS important to make the proper evaluation of o since it has a 
great significance in the outcome of the simplified equations. It is also 
important to realize that the reciprocal of a produces the same values of 
the various parameters given in tables 1 and 1 . as produced by a for a 
given curvature 1/R. 

Figure J shows three diverse situations in which the simplified equa- 
tions can be usefully applied. The locomotive wheel on a rail (fig. 3(a)) 
illustrates an example in which the ellipticity k and radius ratio a are 
less than ). The ball rolling against a flat plate (fig. 3(b)) provides 
pure circular contact (i.e., a » k ■■ 1.00). Figure 3(c) shows how the con- 
tact ellipse IS formed in the ball - outer-ring contact of a ball bearing. 
Here the semimajor axis is normal to the direction of rolling and con- 
sequently a and k are greater than 1. The detailed geometry and the 
values from the calculations that can be made using the simplified formulas 
are given in table t> for each of these configurations. In using these for- 
mulas It IS important to pay attention to the sign of the curvatures. Note 
that the outer race in figure J(c) is a concave surface and so the sign is 
negat i ve. 

SUMMARY OF RESULTS 

An alternative approach has been presented for the classical Hertzian 
solution for local stress and deformation of two elastic bodies in contact. 
^ 111^)11 fled formulas that use curve-fit analysis are given in terms of the 
radius ratio a for the ellipticity K and for the complete elliptic 
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»nte 9 raU, ^ and /, of the first and second kinds, respectively. Thus 
their Inter-dependence can be uncoupled, and solution of the resulting 
transcendental equation avoided. Simplified equations were developed that 
permit a more direct and easy approach to the calculation of the elliptical- 
contact deformation and maximum Hertzian stress. In addition, a curve-fit 
analysis was used to derive a simplified formula for an auxiliary stress 
parameter t as a function of the radius ratio a. This eliminated having 
to solve a cubic equation for t as a function of k. Thus the simplified 
formula for t, together with the simplified formulas for stresses and de- 
formations, permits a direct and easy calculation of the location and magni- 
tude of the maximum subsurface orthogonal shear stress. Therefore the 
elliptical-contact deformation and stresses are presented that are appli- 
cable for any contact ranging from a disk rolling on a plate (radius ratio 
a - 0.03} to a bal 1-on-plate contact (a » 1) to a contact approacning a 
no<n1naI line contact (a * 100) such as a barrel-shaped roller against a 
plate. 
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TABLE S. - PRACTICAL APPLICATIONS FOR DIFFERING CONFORMITIES 


[£• - 2.197x10^ N/ctn^.] 


CONTACT 

PARAMETERS 

wheel on rail 

BALL ON plane 

BALL - OUTER- 
RING CONTACT 

F 

1.00x10^ N 

222.4111 N 

222.4111 N 

■ ax 

50.1900 cm 

0.6350 cm 

0.6350 cm 



0.6350 cm 

0.6350 cm 

^bx 

oe 


-3.8900 cm 

*^by 

30.0000 cm 


-0.6600 cm 

a 

0.5977 

1.0000 

22.0905 

k 

0.7099 

1.0000 

7.3649 

F 

0.7206 

1.0000 

7.1738 

/ 

1.3526 

1.5708 

1.0267 

T 

1.3412 

1.5708 

1.0258 

r 

1.8508 

1.5708 

3.3941 

f 

1.8645 

1.5708 

3.3375 


1.2801 cm 

0.0426 cm 

0.1842 cm 

Oy 

1.2829 cm 

0.0426 cm 
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Figure 1. - Geometry of contacting elastic solids. 
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Figure 2. - Ellipticity parameter and elliptic integrals of first and second kinds as a 
function of curvature ratio. 
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